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ABSTRACT 

We report the detection of a faint optical flash by the 2-m Faulkes Telescope North 
simultaneously with the second of two prompt 7-ray pulses in INTEGRAL gamma- 
ray burst (GRB) 080603A, beginning at trest = 37 s after the onset of the GRB. This 
optical flash appears to be distinct from the subsequent emerging afterglow emission, 
for which we present comprehensive broadband radio to X-ray light curves to 13 days 
post-burst and rigorously test the standard fireball model. The intrinsic extinction 
toward GRB 080603A is high {Av,z — 0.8 mag), and the well-sampled X-ray-to-near- 
infrared spectral energy distribution is interesting in requiring an LMC2 extinction 
profile, in contrast to the majority of GRBs. Comparison of the 7-ray and extinction- 
corrected optical flux densities of the flash rules out an inverse-Compton origin for 
the prompt 7-rays; instead, we suggest that the optical flash could originate from the 
inhomogeneity of the relativistic flow. In this scenario, a large velocity irregularity in 
the flow produces the prompt 7-rays, followed by a milder internal shock at a larger 
radius that would cause the optical flash. Flat 7-ray spectra, roughly F oc u^^-^, are 
observed in many GRBs. If the flat spectrum extends down to the optical band in 
GRB 080603A, the optical flare could be explained as the low-energy tail of the 7-ray 
emission. If this is indeed the case, it provides an important clue to understanding 
the nature of the emission process in the prompt phase of GRBs and highlights the 
importance of deep (i? > 20 mag), rapid follow-up observations capable of detecting 
faint, prompt optical emission. 

Key vifords: gamma-rays: bursts — radiation mechanisms: nonthermal 
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1 INTRODUCTION 

The exact mechanism that produces the prompt radiation of 
a gamma-ray burst (GRB) is still unknown. As a nonthermal 
process, synchrotron and inverse Compton (IC) are the main 
candidates. The former can successfully account for most of 
the afterglow emission evolution, and is naturally expected 
from shock-accelerated electrons. As such, this has also been 
considered for explaining the 7-ray prompt emission itself. 
However, the energy spect rum of a GRB , usually modelled 
with a Band function (,Band et al.lll993l ). is such that the 
typical value for the low-energy photon index (a) violates 
the so-called "synchr otron death line" {a = —2/3) for a siz- 
able fraction of cases l|Preece et al.lll99g : lGuiriec et al.ll2010l : 
iGuidorzi et al.ir201ll ). In addition, the value generally ob- 
served, a « —1, differs from the value of —3/2 expected for 
rapid ly cooling electrons (t he so-called "fast-cooling death 
line" ; (GhiseUini et"al]l2000l) . Although under some assump- 
tions most GRB sp ectra could be reconciled with a syn- 
chrot ron origin (e.g., iLlovd fc Petrosianll2000l : iDaigne et al.l 
[20l2), the question of whether it is the dominant process in 
the GRB production remains unanswered. 

On the other hand, IC has been considered as a pos- 
sible alternative, such as synchrotron self-Compton (SSC; 
iKumar fc McMahonllioOsI ). especially when the p rompt op- 
tical emission is very bright (|Racusin et al.|[2008l ) . IC as the 
source of 7-rays requires a soft component in the infrared 
(IR) through ultraviolet (UV) range for providing the seed 
photons; this in turn means that the second IC component 
peaks in the GeV-T eV range, potentia lly implying an "en- 
ergy crisis" problem (|Piran et al.ll2009l ). Combining prompt 
optical and 7-ray measurements, together with the wealth 
of information derived from the broadband modelling of the 
early-to-late time afterglow, offers a direct way test for IC 
as the mechanism for the GRB production. 

The standard afterglow model (see, e.g.,[Meszaros"2006' 
for a review) is rather successful in modelling the temporal 
and spectral evolution of GRB afterglows. However, models 
often require modifications, such as energy injections from 
long-lived internal engines or density enhancements in the 
surrounding medium, and for a sizable fraction of cases, even 
these options cannot provi de a fully satisfactory explanation 
(e.g., iMelandri et"al](2008l ). 

Long-duration GRBs are also probes of the interstel- 
lar and intergalactic medium, and of cosmic star forma- 
tion history up to redshift z ~ 8 (ISalvaterra et al] l2009l : 
iTanvir et al. 20091). potentially exploring the re ionization 
epoch l|Kistler et al.1 [20091 : [Robertson et al.ll2010l ). Spectral 
energy distribution (SED) studies and spectroscopic obser- 
vations can help to shed light, for instance, on the redshift 
evolution of dust, gas content, and metallicity of the host 
galaxies o f GRBs as well as of t he local region within the 
host (e.g.. [Prochaska et al.|[2007l ). In addition, they help to 
identify the crucial physical parameters which favour the 
production of GRBs. Dust-extinction modelling for a sample 
of GRB afterglows with well-sampled SEDs has shown that 
most cases can be describe d with Small Magellanic Cloud 
(SMC) profiles (|Kann et al.i 2010) having little evidence for 
the 2175 A bump seen i n the Milky Way extinction curve, ex- 
cept f or very few cases jKriihler et al.ll2008l: Prochaska et alj 



I2OO9I : lElfasdottir et all 120091 : iPerlev et al.ll201ll )~ 

No direct link has been found between the properties of 



the prompt emission and those of the circumburst environ- 
ment surrou nding the GRBs and of the host galaxy, such as 
metallicity (jLevesgue et al.ll2010l ). Nevertheless, a detailed 
picture of the properties that can be derived from the broad- 
band afterglow modelling — the dust content and features 
along the sightline to the GRB within the host galaxy — 
are crucial to provide a self-consistent description of the en- 
tire GRB phenomenon, and for unveiling the yet unknown 
connections between the GRB itself and its birthplace and 
(to some extent) progenitor. 

This paper reports comprehensive analysis and dis- 
cussion of the multi-wavelength dataset collected on the 
long-duration GRB 080603A detected by INTEGRAL 
jWinkler et al.l 12003 ) in light of the current standard fire- 
ball model. This GRB provides an ideal test bed because 
it had an optical flash simultaneous with the prompt emis- 
sion, and we recorded the broadband afterglow SED and 
its evolution. Our dataset includes INTEGRAL data of the 
7-ray prompt emission itself, as well as multi-filter photo- 
metric and spectroscopic data of the near-infrared/optical 
afterglow and of the host galaxy. In addition, we analysed 
the X-ray afterglo w data, discovered in th e 0.3-10 keV band 
with Smft/XRT (|Sbarufatti et alj|2008bl) . from 3 hours to 
7 days after the burst. We also include data from the Very 
Large Array (VLA), taken from 2 to 13 days post burst, in 
which the radio afterglow was detected. 

Throughout the paper, times are UT and are given 
relative to the GRB onset time as observed with INTE- 
GRAL, which corresponds to June 3, 2008, 11:18:11 UT. 
The convention F(v,t) oc t'" is followed, where the 
energy index j3 is related to the photon index by F = 
(3-1-1. We adopted the standard cosmological model: Hq = 
71kms-^Mpc-\ Qa = 0.73, D.m = 0.27 (jSpergel et all 
I2OO3I ). 

All of the quoted errors are given at 90% confidence 
level for one interesting parameter (Ax^ ~ 2.706), unless 
stated otherwise. 



2 OBSERVATIONS 

GRB 080603A was detected with the INTEGRAL/IBIS in- 
strument and localised in real time by the IN TEGRAL 
Burst Alert System (IBAS; iMereghetti et aLll2003i ) with an 
accuracy of 3.2'. The 7-ray prompt emission in the 20- 
200 keV energy band lasted about 150 s. A quick-look 
analysis gave a peak flux of 0.5 ph cm~^ s~^, a fluence 
of about 10~^ erg cm~^, and burst coordinates a(J2000) 
= 18'' 37"'38!2, J(J2000) = -h62°44'06" with an error radius 
of 2' (|Paizis et al.ll2008l ). 

The Faulkes Telescope North (FTN) promptly reacted 
to the IBAS alert and began observing at 11:19:51, which 
was 100 s after the burst onset time and while prompt 7-ray 
emission was still ongoing. The automatic GRB pipeline did 
not identify any candidate; consequently, a pre-programmed 
BVRi' observation sequence wi th increasingly longe r expo- 
sure durations was carried out (|Guidorzi et al.ll2006l ). How- 
ever, visual inspection of the frames revealed the presence of 
an uncatalogued, variable object at a(J2000) = 18''37'"38'!1, 
5(J2000) = -F62°44'39'.'4 with R = 19.6 mag at t = 7.37 min 
calib rated against nearby USNOB-1.0 stars (|Gomboc et al.l 
l2008h . FTN observations continued for 3 hours. 
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The GRB alert syste m of the Katz man Automatic 
Imaging Telescope (KAIT; iLi et al.l |2003| ) at Lick Obser- 
vatory also promptly reacted to the INTEGRAL alert and 
independently detected the optical counterpart at a position 
consistent with that of the FTN, reporting I = 18.7 mag at 
10.7 min in unfiltered and 7-band images (jChornock et alJ 
|2008| ). Other rob otic telescopes also rep orted the discovery 
of the afterglow (|Milne fc Updike 2008^. 

We began spectroscopic observations of the opti- 
cal afterglow with the GMOS dual spectrometer at the 
Gemini-North 8-m telescope starting at time 13:24, iden- 
tifying several absorption features at a common redshift 
oi z = 1.68742, in agreement with preliminary reports 
l|Perlev et al.ll2008ah . 

We kept monitoring the evolution of the optical after- 
glow with the 1.34 m Schmidt telescope of the Thiiringer 
Landessternwarte (TLS), the Liverpool Telescope (LT), 
the Keck-I 10-m telescope, the Zeiss-1000 and AZT- 
11 telescopes at the Crimean Astrophysical Observatory 
(CrAO), and the AZT-33IK telesco pe at Sayan Ob- 
servatory up to 4 days post b urs t (iKann et al.' 200^; 
Rumvantsev fc Pozaiienkol l2008allbl: IPerlev et al. 2008bi; 



600 



Klunko fc PozanenkolbOOanRumvantsev et aLlboOSl ). From 



19 to 22 hours post burst we observed the afterglow at near- 
infrared (NIR) wavelengths with the Peters Automated In- 
frare d Imaging T elescope (PAIRITEL) through JHK fil- 
ters l|Miller et al.l 120081. Our last detection of the optical 
afterglow was obtain ed with th e Low Resolution Imaging 
Spectrometer (LRIS; lOke et al] Il995t ) 4 days post burst. 
Two months after the burst we observed the same field 
with Keck/LRIS using u'g'RI filters and detected the host 
galaxy; these observations were subsequently used to sub- 
tract the host contribution when the afterglow became com- 
parably faint. 

Swift/XRT began observing GRB 080603A from 2.9 
to 5.9 hours and reobserved it from 2.5 to 7.0 days af- 
ter the burst. The X-ray afterglow was clearly identi- 
fied at position a(J2000) = 18''37™38"06j_5{J2000}_= 
+62^4' 40'.' 1 with an error radius of 1'.'9 (|Sbarufatti et all 
l2008al ). During the same time interv al Swift/UVOT detected 
the fading afterglow in the V band (|Kuin fc Manganoll2008l : 
ISbarufatti et al.ll2008bl ). 

Finally, we discovered the radio counterpart with the 
VLA at 4.86 and 8.46 GHz, initially at the latter frequency 
with a possibl e detection at 1.9 days at a flux density of 
116 ± 41 /xjy llChandra fc FraillbOOSal). The det ection was 
confirmed at 3.95 davs ijChandra fc Frailll2008bh . Observa- 
tions continued as late as 13 days post burst. 

The Galactic reddeni ng along the d i rectio n to the GRB 
is Eb-v = 0.044 mag ()Schlegel et al.lll998l l. The Galac- 
tic extinction in each filter has been estimated through the 
NASA /IPAC Extragalactic Database extinction calculatoiQ. 
Specifically, the extincti on in each filt e r is de rived through 
the parametrisation by ICardelli et al.l (Il989l ): A^i = 0.23, 
Ab = 0.19, Ag, = 0.18, Av = 0.14, A^' = Ar = 0.12, 
A,, = 0.09, Ai = 0.08, Aj = 0.04, Ah = 0.025, and 
Ak = 0.02 mag. 



http: / /nedwww. ipac.caltech.edu/forms/calculator. html. 
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Figure 1. INTEGRAL light curve in the 20-200 keV energy band 
(left-hand ordinate axis). The shaded areas display the time in- 
tervals of the first optical measurements with the FTN; the corre- 
sponding ii-band flux densities are shown with filled circles (right- 
hand ordinate axis) . The dashed line shows the best- fitting models 
of the 7-ray pulses as obtained with the model by iNorris et al.l 
l|2005l l. 



3 DATA REDUCTION AND ANALYSIS 
3.1 Gamma-ray data 

Figure [1] shows the 20-200 keV background-subtracted time 
profile of GRB 080603A recorded by the IBIS/ISGRI de- 
tector (|Lebrun et al.ll2003l ). The profile consists of two very 
similar pulses of duration 30-s, peaking at 14 s and 114 s. 
A combination of two fast- rise-exponen t ial-de cline (FRED)- 
like pulses as modelled bv lNorris et al.l (|2005l ) gives a satis- 
factory result {x^/doi — 94.4/81), as shown by the dashed 
line in Figure [T] The parameters used are the peak time 
tpeak, the peak intensity A, the rise and decay times Tr and 
Td, the pulse width w, and the asymmetry k. Their best- 
fitting values are reported in Table [T] Apart from the peak 
of the second pulse, which is roughly twice as intense as that 
of the first, the two pulses share very similar temporal prop- 
erties: rise and decay times around 7 and 20 s, respectively, 
with a corresponding decay-to-rise ratio around a fac tor of 
3, very typical of classical FREDs (|Norris et al.l ll996''). 

Two time-integrated spectra, one for each pulse, show 
no evidence for spectral evolution: a simple power law can 
fit both spectra with a 7-ray photon index F^. « 1.6. Table[2] 
reports the best-fitting spectral parameters. The 20-200 keV 
total fluence is (1.1 ± 0.2) x 10~^ erg cm~^ , in agreement 
with preliminary reports jPaizis et al.ll200^ . The value of 
r^, = 1.6 lying between the typical low-energy and high- 
energy photon indices of GRB prompt emission spectra (e.g., 
Kaneko et al.,.2006 ; Sakamoto ct al. 2009) suggests that the 
peak energy, Ep, is likely to lie within the 20-200 keV en- 
ergy band. In the cont ext of the sample of GRBs detected 
with ISGRI and BAT l|Vianello et al.fl2009l ). the fiuence of 
GRB 080603A makes it a typical burst. 

The corresponding flux-density curve shown in Figure S] 
was found to refer to 84 keV; this is the energy at which the 
energy spectrum with Z^^. = F^. — 1 = 0.6 has the same value 
as that averaged over the 20-200 keV range. 

Despite the unknown value of Ep , we can provide a con- 
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Table 1. Best- fitting parameters of the time profile of tiie 7-ray pulses as seen in the 20-200 keV band. 



Pulse 


^peak 


A(at 84 keV) 






w 


k 




(s) 


(/^Jy) 


(s) 


(s) 






1 


13.7 ± 1.9 


18 ±2 


8.0 ±2.0 


21.8 ± 3.5 


29.8 ±3.8 


0.46 ±0.12 


2 


113.8 ±0.9 


41 ± 2 


6.2±1.1 


19.4 ± 1.6 


25.7 ± 1.8 


0.51 ±0.08 



Table 2. Best-fitting parameters of the energy spectra of the 7-ray prompt emission in the 20—200 keV band. The model is a power-law 
and Fj is the photon index. 



Pulse 


Time interval 
(s) 




Average flux 
(10~* erg cm~^ s^-"-) 


xVdof 


Fluence 
(10"'^ erg cm~^) 


1 
2 


3-38 
100-140 


1.6 ±0.2 


1 4+0'14 
^■^-0.6 
1 n+0.16 
^•^-0.6 


60/64 
39/30 


4.1 ± 1.3 
6.7± 1.5 



servative estimate of the isotropic-equivalent radiated en- 
ergy Eiso in the GRB rest-frame 1-10* keV energy band: 
we assume that Ep lies either within or close to the 20- 
200 keV energy range. In the former case, we use the loga- 
rithmic average, Ep = 60 keV, while in the latter case we 
consider the values 10 keV and 400 keV as the lower and 
upper boundary, respectively. These values correspond to a 
0.3 logarithmic shift from the corresponding boundary, the 
logarithmic bandwidth being 1. In calculating the fluence in 
the rest-frame 1-10* keV band, the K-correction factor is 
2.8 ± 0.8, where the uncertainty accounts for the different 
Ep assumed and where we adopte d the typical Band f unc- 
tion with ob = -1 and I3b = -2.3 (|Kaneko et al.ll2006l '). As 
a result, we estimate Eiso = (2.2 ± 0.8) x lO^'^ erg and the 
intrinsic peak energy Ep^i — I6OI130 keV; the se values are 
broadly c onsistent with the Ep^i-Eiso relation (|Amati et al.l 
|2002| : lAmati 2010), although the poor accuracy on Ep^i is 
not very constraining. 



3.2 X-ray data 

The Swift/XRT began observing GRB 080603A on 2008 
June 03 at 14:11:19, about 10.4 ks after the burst, and ended 
on 2008 June 10 at 11:44:56, with a total net exposure of 
17.8 ks in photon counting (PC) mode spread over 6.9 days. 
The XRT data were processed using the FTOOLS software 
package (v. 6.7) distributed within HEASOFT. We ran the 
task XRTPIPELINE (v. 0.12.1), applying calibration and stan- 
dard filtering and screening criteria. Data were acquired only 
in PC mode due to the faintness of the source. Events with 
grades 0-12 were selected. The XRT analysis was performed 
in the 0.3-10 keV energy band. 

Source photons were extracted from a c ircular region 
centred on the final XRT position (^Sbarafatti et al. 2008]a) 
and with a radius of 20 pixels (1 pixel = 2'.' 36), and were 
point-spread function (PSF) renormalised. Background pho- 
tons were extracted from nearby circular regions with a total 
area of 22.7 x 10^ pixels away from any source present in the 
field. No pile-up correction was required because of the low 
count rate (< 0.1 count s^^) of the source from the begin- 
ning of the XRT observations. When the count rate dropped 
below ^ 10~^ count s~^, we made use of ximage with the 
tool SOSTA, which corrects for vignetting, exposure varia- 



tions, and PSF losses within an optimised box, using the 
same background region. 

We extracted the 0.3-10 keV energy spectrum in the 
time interval from 10.4 to 21.1 ks; later observations did not 
allow us to collect enough photons to ensure the extraction of 
another meaningful spectrum. Source and background spec- 
tra were extracted from the same regions as those used for 
the light curve. Spectral channels were grouped so as to have 
at least 20 counts per bin. The ancillary response files were 
generated using the task xrtmkarf. Spectral fitting was 
performed with xspec (v. 12.5). The spectrum can be mod- 
elled with an absorbed power law with the combination of 
XSPEC models wabs zwabs pow, based on the p hotoelectric 
cross section by iMorrison fc McCammonI (119831 ). Results of 
the best-fit parameters are reported in Table|6l The Galactic 
neutral hydrogen column density along the GRB direction 
was fixed to the value determi ned from 21 cm line radio sur- 
veys: N§f = 4.7x 10^° cm-2 (|Kalberla et allbOOSl ). The ad- 
ditional X-ray absorption, modelled in the GRB rest frame, 
was found to be A^hi.z = 6.6l4 g x 10'^^ cm~^, very typ i- 
cal of X-ray afterglow spectra (e.g.. [Campana et akllioiol ). 
The X-ray photon index in the 0.3-10 keV energy band is 
Fx =2.3 ±0.3. 

The X-ray unabsorbed flux light curve was derived from 
the rate curve by assuming the same counts-to-energy fac- 
tor (5.4 X 10~^^ erg cm"'^ s~^ count"^) obtained from the 
spectrum described above. This implicitly relies on the lack 
of strong spectral evolution from ^ 10 ks onward; although 
such an assumption cannot be proven due to the paucity 
of photons at late times, this i s in agreement wit h what is 
observed for most GRBs (e.g.. lEvans et al . 2009). Finally, 
the flux-density curve shown in Figure |4] was calculated at 
2.4 keV, the energy at which the energy spectrum with 
Px ~ Fx — 1 = 1.3 has the same value as that averaged 
over the 0.3-10 keV range. 

3.3 Infrared/optical data 

The FTN carried out robotically triggered observations be- 
tween 100 s and 190 min. During the detection mode, con- 
sisting of the flrst 3 x 10 s frames in the R band, the op- 
tical afterglow was too faint to b e automatically iden tified, 
so the GRB pipehne LT-TRAP (jCuidorzi et al.ll2006l ) trig- 
gered the multi-filter {BVRi') observation sequence with 
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increasingly longer exposures. However, a quick visual in- 
spection of the data led to the identification of an uncata- 
logued and variable sou rce proposed to be the afterglow can- 
didate jGomboc et al.ll2008 ,y Our best estimate for the opti- 
cal afterglow position is a(J2000) = 18''37™38=05, J(J2000) 
= -|-62°44'39'.'4 with an error radius of 0'.'5, and it lies within 
the final XRT error circle. 

The afterglow observations with the FTN occurred dur- 
ing the onset of the second 7-ray pulse; concurrently, a faint 
optical flash with an R magnitude varying within the 20-21 
range was observed (Fig. The flash was soon followed by 
a steep rise and a broad plateau around 10^ s, at the end 
of which a smooth transition to a typical power-law decay 
with index around 1 took place (Fig. 3)). 

Later observations were carried out with the LT from 
15.5 to 17.1 hours with the Sloan Digital Sky Survey (SDSS) 
r' and i' filters, as well as with the FTN from 21.4 to 
26.4 hours in i' . 

Calibration of the BVRr'i' frames was performed 
by comparing with the magnitudes of four nonsaturated 
field stars. The corresponding zero points were determined 
throu gh the observations of iLandolt (199^) field stars for 
which Smith et al.l l|2002l 'l provide an SDSS calibration. In 
both cases the zero points were stable during the night, 
showing fluctuations as large as 0.02-0.03 mag in the worst 
cases. Finally, we corrected for the airmass. Both aperture 
and PSF photometry was systematically carried out using 
the Starlink GAIA softwarqj, making sure that both gave 
consistent results within the uncertainties. Magnitudes were 
conve r ted in to flux densities (/iJy) following iFukugita et al.l 
i|i995l . Iigge'v Results are reported in Table [S] magnitudes 
are corrected for airmass, while flux densities are also cor- 
rected for Galactic reddening. 

KAIT observations began a,t t — 253 s in the V, I, and 
unfiltered bands; the pre-programmed exposure durations 
generally increased with time. The earliest firm unfiltered 
detection of the optic al afterglow occurred at t ~ 11 min 
l|Chornock et al.ll2008l ). independent of the FTN detection. 
Successive stacked frames also gave a later detection with 
the / filter, while the V filter only provided either upper 
limits or very marginal detections, as reported in Table [H 
For the V filter, w e also considered t he UV OT photometric 
points provided bv lSbarufatti et al] (|2008bl ). We used some 
of the four FTN field stars to cahbrate the KAIT field; for 
the unfiltered frames we adopted the zero point of the R 
band. The last useful frame obtained with KAIT was ac- 
quired at 20 min. Despite the large uncertainties, both the 
measured values and upper limits are in agreement with the 
contemporaneous values obtained with the FTN, as shown 
by Figure |4l 

At a midpoint time of 1.55 days we observed 
GRB 080603A with the 1.34-m Schmid t TLS, obtaining a 
total of 6 X 600 s images in the R band l|Kann et al.ll2008l ). 
Stacking all six frames, the afterglow is detected. In order 
to subtract the contribution from the extended object later 
identified as the host galaxy and properly account for the 
crowded field, another frame of the same field was taken on 
2008 August 26 and used for image subtraction with the ISIS 



packag^l l|Alard fc Luptoiilll998l ). The brightness, estimated 
with both aperture ph otometry and SExtractor (v. 2.5.0; 
iBertin fc Arnoutslll99d ). is 7? = 22.1 ± 0.3 mag. 

Other late-time observations with R filters were 
obtained with the 1.25-m AZT-11 (at 0.36 days) and 
1-m Zeiss-1000 (ZIOOO, at 1.6 and 3.6 days) telescopes 
of the Crimean Astrop h ysical Obs ervatory (CRAQ; 
^Rur nvantsev fc Pozanenkq l2008al lbl: iRumvantsev et all 
a004 ). and with the 1.5-m AZT -33IK telescope (at 
0.30 days) of the Sayan Observatory (|Klunko fc Pozanenkd 
I2OO8I ). In all cases the afterglow was detected, except for the 
last observation at 3.6 days, which provided an upper limit 
of _R > 22.9 mag. For the same reasons as in the case of the 
TLS frame, we had to correct the measured R magnitude 
at 1.6 days for the host-galaxy contribution. Given the lack 
of late-time images, we merely subtracted the host-galaxy 
fiux contribution as estimated with Keck (see below); this 
turned into a shift of 0.2 mag in R, comparable with the 
uncertainty affecting the measurement itself, which in the 
end was R = 22.52 ± 0.25 mag. 

The 1.3-m PAIRITEL started observing the afterglow 
of GRB 080603A at 18.9 hours with JHK^ filters. Photo- 
metric calibration was done against seven nearby 2MASS 
stars; magnitudes were estimated with both aperture and 
PSF photometry under GAIA. The NIR afterglow counter- 
part is clearly detected in all filters in two mosaic frames 
centred at 0.82 and 0.95 days with 2822 s and 4363 s total 
exposures, respectively. Our estimat es agree within un cer- 
tainties with the preliminary results (iMiller et al.ll20()8l ). 

3.3.1 Late-time host galaxy observations 

We used the Keck LRIS to observe GRB 080603A at two 
different epochs. The first run was taken on 2008 June 7, 
between times 12:31 and 12:47 with R (total exposure 690 s) 
and g' (total exposure 785 s) filters. The average airmass 
was 1.38 and the D560 dichroic was used. The afterglow 
was clearly detected with both filters at the FTN position. 

The same field was reobserved on 2008 August 2, with R 
(total exposure 930 s) and g' (total exposure 1110 s) filters. 
The average airmass was 1.58 and the D560 dichroic was 
used. The next night, we observed with the u' and I filters 
(total exposure 720 s each), using the D680 dichroic, with 
an average airmass of 2.00. 

Calibration was based on the identification of seven 
common field stars having R = 20-23 mag, bright enough 
to be accurately measured in the stacked images of the 
FTN-fLT and faint enough to avoid saturation in the Keck 
images. Transformations between the SDSS and Johnson- 
Co usins systems for t his set of faint stars were done follow- 
ing (jOTdTer^l] (|2006l ). The scatter in the zero point of each 
filter was incorporated into the uncertainties. 

In order to subtract the contribution of the host galaxy 
from the afterglow flux observed at 4.1 days, we performed 
image subtraction with ISIS. Figure [5] shows the result in 
the R band: the afterglow plus host at 4.1 days, the host 
two months later, and the difference between the two are 
shown in the left, middle, and right panels, respectively. 
From image subtraction, we estimated R — 24. 17 ±0.14 and 



^ http:/ /starlink.jach. hawaii.edu/starlink. 



^ http:/ /www2. iap.fr/users/alard/package. html. 
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Figure 2. Keck-LRIS /?-band images of tiie crowded field of GRB 080603A. The circle is our best position of tfie optical afterglow 
obtained from early FTN frames with 0''5 error radius. Left panel: taken at 4.1 days; the afterglow is still detected. Middle panel: taken 
two months after the burst; the host galaxy is clearly visible within the error circle, while the afterglow is no longer detectable. Right 
panel: subtracted image. 



Table 3. Host-galaxy magnitudes. 



Fihcr 


Aoff 

(A) 


Magnitude 


Corrected^ 


u' 


3450 


25.6 ± 0.4 


25.4 ±0.4 


g' 


4731 


24.78 ±0.16 


24.60 ±0.16 


R 


6417 


24.03 ±0.12 


23.91 ±0.12 


I 


7599 


23.98 ±0.15'' 


23.89 ±0.15*^ 



^ For Galactic extinction. 

^ Calibration was done against i' magnitudes of the field stars, 
corresponding to Aoff = 7439 A. 

-Rhost = 24.03 ± 0.12 mag for the afterglow and the host, re- 
spectively. Correspondingly, we found g' — 24.95 ± 0.20 and 
Shost = 24.78 ± 0.16 mag. We also used the late-time R and 
I frames to subtract the host contribution from the obser- 
vations performed with the FTN and LT at t ~ 1 day, when 
the afterglow brightness was around 21 mag in the same fil- 
ter. In the two cases this turned into a shift of ~ 0.1 mag, 
similar to the corresponding statistical uncertainties. For the 
same reasons, we did not correct the comparably faint after- 
glow magnitudes of the optical flash seen in the first FTN 
images, because of their relatively large uncertainties. 

The u'-band magnitude of the host was calibrated with 
observat i ons o f three standard stars in the PG 0231+051 
iLandoltl (|l992l ) field taken the same night with an airmass 
of 1.0. Assuming an extinction coefficient in the range 0.36- 
0.4, the zero point for the u'-band filter was estimated to be 
Z — 27.8 ± 0.3 mag. Table [3] reports the photometry of the 
host galaxy. 

The host centroid, as determined with SExtractor, is 
a(J2000) = 18''37'"38=03, 5(J2000) = +62°44'39'.'0, which 
is 0'.'4 away from the afterglow position. This angular off- 
set corresponds to a projected distance of 3.4 kpc. Taking 
into account the uncertainty in the afterglow position, an 
upper limit of 6 kpc is more conservative, in agreement 
with the typical p rojected offsets of long-duration GRBs 
l|Bloom et ai.ll2002l ). 



reduced with the LowRedux pipelineQ developed by J.X.P. 
and J. Hennawi. This custom software bias subtracts, flat- 
fields, and optimally extracts the spectra. The two exposures 
were coadded, weighting by signal-to- noise ratio (S/N); the 
resultant spectrum has a S/N exceeding 80 per pixel over 
the range A ^ 5000-8000 A. 

The spectrum reveals a series of very strong rest-frame 
UV transitions at ^ = 1.68742 that includes Al III AA1854, 
1862, Zn II A2026, Mg II, and a plethora of Fe II res- 
onance and fine-structure transitions (Fig. [3)|. The detec- 
tion of the fine-structure transitions uniquely identifies this 
gas (and redshift) with the host galaxy of GRB 080603 A 
(|Prochaska et al]|2006l ). We have also detected two inter- 
vening Mg II systems along the sightline at z = 1.2714 
and z = 1.5636. The former marks yet another example 
of a strong (equivalent width W2796 > 2 A) Mg II system 
along a GRB sightline JProchter et aLlbOod l. Table [1 lists 
the equivalent-width measurements (from Gaussian fits to 
the absorption lines) for all of the features detected in our 
spectrum. Note that the total equivalent width is reported 
for transitions that are severely blended. 



3.4 Radio data 

GRB 080603A was observed with the VLAQat four epochs in 
the 8.5 GHz band and at 2 epochs in the 4.9 GHz band. Our 
observations spanned from 2008 June 5 until June 16. That 
on June 5 was made in the VLA C-configuration, whereas 
the later observations were made in the DnC-configuration. 
We adopted the VLA calibrator J1835+613 for phase cali- 
bration at both frequency bands. 

The data were analysed using standard data-reduction 
routines of the Astronomical Image Processing System 
(AIPS). 3C 286 was used for the fiux calibration. We had 
a possible detection (~ 2.8cr) in the first observation at 

8.5 GHz; it was confirmed in the subsequent observations. 
The results are reported in Table H) 



3.3.2 Spectroscopy 

We initiated observations of the GRB afterglow with the 
Gemini/GMOS dual spectrometer starting at time 13:24 for 
a series of two 1200 s exposures using the R400 grating, off- 
set by ~ 5 A to fill in the detector gaps. The data were 



* http: / /www. ucolick.org/ ~xavier/LowRedux/index.html. 
^ The NRAO is a facility of the National Science Foundation 
(NSF), operated under cooperative agreement by Associated Uni- 
versities, Inc. 
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Figure 3. Gemini/GMOS spectrum of the afterglow taken at 2.1 hours post burst. 
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Figure 4. Panchromatic light curves of the prompt and afterglow emissions. 
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Table 4. 
VLA. 



Flux densities of the radio afterglow obtained with the 



Table 5. Multi- filter light curve best-fitting parameters. 



UT Date 



At 

(days) 



(GHz) 



(A'Jy) 



June 05.39 


1.92 


8.46 


116 ±41 


June 07.42 


3.95 


8.46 


154 ± 28 


June 08.24 


4.77 


4.86 


112 ±45 


June 08.26 


4.79 


8.46 


230 ± 29 


June 16. f 9 


12.72 


4.86 


186 ± 49 


June 16.21 


12.74 


8.46 


70 ±42 



4 RESULTS 

4.1 Multi-band light curves 

Figure |4] shows the 7-ray prompt and broadband afterglow 
light curves. In modelling the data we initially allowed colour 
change between 5 x 10^ and 10^ s, i.e. the best-sampled in- 
terval, by fitting the data with the sum of a fast-decaying 
component and a slow-decay component. The colour change 
between the two was 0.225 ± 0.344, so less than 1-a. De- 
pending on whether or not we allow colour change, the re- 
sults do not change within uncertainties. Given the apparent 
lack of chromatic changes, to better constrain the evolu- 
tion we simultaneously fitted the light curves of the various 
bands (except at the radio wavelengths) with the same func- 
tion, only allowing different normalisations and no colour 
change. We modelled the different power-law regimes with 
the smoothly broken p ower-law model parametrisation by 
iBeuermann et aU l|l999l l. Furthermore, given the clear pres- 
ence of a break in the R/r' and B/g' curves around 10^ s, 
we allowed a further achromatic break. The final model is 
described by Equation ([l}. 

{a2-Q3)/n2 



Fit) = Fo 



1+ 



+ 



l/n 



(1) 



The best-fitting solution was found through minimisation of 
the overall x^, resulting from the sum of the total values 
of the individual light curves with respect to each corre- 
sponding model. The free parameters are the following: cti, 
02, and Qf3 are the power-law indices during the initial rise, 
the following decay, and the final (> 10'' s) decay, respec- 
tively. The two break times are thi and th2, while n and 712 
are the smoothness parameters regarding the first and sec- 
ond breaks, respectively. Only 722 could not be determined 
from the fit because of the sparseness of the data around the 
final break and it was therefore fixed to 10, so as to give a 
rather sharp break. The normalisation term is represented 
with Fo in Equation ^ , although in practice the free param- 
eter we used for each profile was the fiux density calculated 
at a fixed reference time (we chose 1.5 x 10* s, close to the 
time of most observations). The purpose of this choice is to 
limit the effects of the strong correlations between some pa- 
rameters, such as ai and tbi, on the determination of the 
confidence intervals for each parameter. 

Concerning the late-time g' point, we converted its fiux 
density into that corresponding to the close-in-frequency B 



band by the factor / = {vB/vg 



0.9, where Pa 



Parameter 


Value 


Unit 


ai 


-3.6 ± 1.7 




*bl 


610t^™ 


s 


«2 


0.99 ± 0.07 




n 


30+?'^? 




UK 


505+1^3 


^Jy 


riH 




MJy 


nj 




MJy 


rii 


73 ±3 


MJy 


rir 


42.0 ± 1.5 


MJy 


nv 


OQ Q-t-6.0 

•^^•^-5.0 


MJy 


riB 


18.0 ±0.9 


MJy 


n-x 


0.22 ±0.03 


^Jy 


th2 




10^ s 


013 


1 7+0.4 




xVdof 


102.4/86 





visible wavelengths (Section I4.2|l . This allows us to better 
constrain the late-time break, in addition to that offered by 
the J?-band curve. The same correction was applied to the 
few r' flux-density points to shift them into the R band. 

This correction introduces a negligible systematic un- 
certainty. We verified this by alternatively fitting a SED 
derived by treating all the filters separately, and the best- 
fitting parameters and corresponding uncertainties turned 
out to fully agree with the values obtained by merging the 
filte rs abov e mentioned and presented in Section 14.21 (see 
also iKann et al.l I2OI1I ) . 

The best-fitting models are shown in Figure |3] and the 
corresponding best-fitting parameters are reported in Ta- 
ble O where the normalisation terms are expressed as flux 
densities calculated at 1.5 x 10* s. 

The peak time, tp, is a function of the free parameters 
as expressed by Equation ((2]): 

l/n{a2—ai) 

tp = thl 



i—y 

V 02 / 



1575tg° 



(2) 



is the observed spectral index of the SED at the observed 



The uncertainty in tp was calculated through error propa- 
gation, taking into account the covariance of parameters. 
Equation ([2]) is exact only when no further breaks are 
present — that is, when tb2 ~ 00 (or, equivalently, 03 — 0:2). 
In practice, it still holds provided that fb2 S> tp, as in this 
case. 

From the flt we excluded the earliest points at t < 155 s, 
connected with the optical flash observed contemporane- 
ously with the last 7-ray pulse (see Fig. [4|. We accounted 
for the presence of some degree of variability around the 
best-sampled curves, particularly around the broad peak at 
tp, by adding in quadrature a systematic error to the sta- 
tistical ones (7%, 4%, and 8% for the B, V, and i' curves, 
respectively), so as to yield a satisfactory goodness of the 
fit: xVdof = 102.4/86 (p-value of 11%). The reason for the 
additional errors was to avoid the risk of underestimating 
the fit parameters' uncertainties. We also note that the na- 
ture of this additional systematic scatter cannot be entirely 
ascribed to unaccounted variability of the zero points; at 
least a few percent must genuinely characterise the afterglow 
light curve. Indeed, the largest deviations from the models 
are seen to occur simultaneously and to correlate in all of 
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the more densely sampled curves: B, R, and i' . Such devi- 
ations are mostly observed within the interval 500-1000 s 
(final part of the rise) and around 10* s, at the beginning of 
the decay. 

We alternatively adopted a different log-likelihood func- 
tion, one more general than that connected with the of 
equation ((l| in that it treats the additional systematic errors 
as free parameters. In practice, this approach does not pro- 
vide any noticeable difference in the best-fitting parameters 
and uncertainties, and the physical implications discussed 
in Section [5] are completely unaffected. The only difference 
concerns the systematic error affecting the B-filter values, 
for which the alternative log-likelihood provides a systematic 
error compatible with zero within uncertainties. In any case, 
this does not affect the best-fitting model to any noticeable 
degree. 

Admittedly, because of the paucity of NIR points (only 
two, very close in time in each filter), the fitting mod- 
els in Figure |4] at the corresponding wavelengths assume 
an achromatic evolution extended to the NIR bands. At 
first glance this may seem too arbitrary, as there have 
been GRBs, particularly those with SEDS that are tem- 
porally well resolved, for which some chro matic evolu- 
tion w as observed, such as GRB 61126 (jPerlev et al.l 
l2008d). GRB 07102 5 (iPerlev et al.ll2010^ and GRB 080319B 



I Bloom et ahl l2009l : [Racusin et all 120081 ). Nonetheless 



note that no chromatic evolution was required by compa- 
rably early-time obs e rvations of several other GRBs (e.g., 



Kriihler et al.1 l2009l: . . 
2010l : ICovino et al.ll201 



Nvse wander et al. . 
Ol: IPerlev et al.ll201ll i 



20091 : lYuan et alj 



Such an assumption has also been made implicitly for 
the X-ray data; these are too sparse to be modelled inde- 
pendently. Our results show that the X-ray light curve can 
be described with the same rigid model as the NIR/optical, 
but of course other possibilities cannot be ruled out. Fur- 
thermore, the two latest X-ray points in Figure |4] assume 
the same spectral shape as that observed around 1.5 x 10* 
due to paucity of X-ray photons observed after 10^ s (Sec- 
tion [3I2J. Although the X-ray photon index of 2.3 observed 
at 10* s is not expected to significantly evolve, the assump- 
tion of no X-ray spectral evolution from 10* s to 10^ s must 
be pointed out. 

The radio data clearly show a different behaviour: the 
peak of the afterglow spectrum crosses the radio bands a 
few days later. Because of this, they were not considered in 
the achromatic modelling, but deserve a dedicated analysis 
in the framework of the standard afterglow model discussed 
in Section [STTI 



4.2 Spectral energy distribution 

We derived an SED from the multi-filter light-curve fitting. 
Although this is based on the achromatic evolution and, as 
such, does not refer to a particular epoch, we considered the 
reference time Uei = 1.5 x 10 s as the most representative 
of it: this is the midpoint time of the X-ray data, when the 
high-energy end of the SED does not rely on any assumption. 
Figure [5] displays the resulting GRB rest-frame SED. 

We considered two different models, either a simple 
power-law and a broken power law with a cooling break, 
A/3 — 0.5, combined with three di fferent exti nction profiles 
according to the parametrisation of|Pel (|l992l ): SMC, LMC2 




Figure 5. Rest-frame SED at 1.51 X 10'' s. The solid line shows 
the best-fitting power law with a LMC2 dust-extinetion profile. 
The dashed line shows the same unextinguished power law. 



(for the Large MageUanic Cloud), and Milky Way (MW). 
Note that while Pel's measurements of the average MW and 
SMC extinctio n curves are general ly consistent with recent 
estimates (e.g.. [Gordon et "al]l2003l ), the LMC implementa- 
tion is representative only of the area around 30 Doradus 
(the LMC2 supershell); hence, we denote this curve as LMC2 
in this work. The average LMC curve is much more similar 
to that of the Milky Way. 

Table[6]reports all of the results of fitting the SED with 
different extinction profiles and/or different models at differ- 
ent energy ranges. The uncertainties on the best-fitting pa- 
rameters include the dependence of the effective frequency 
of each filter on the folded model. Interestingly, only the 
LMC2 profile can satisfactorily account for the NIR/optical 
SED, possibly including the presence of the 2175 A bump. 
The dust content is remarkable: the rest-frame extinction 
is Av,z ~ 0.80 ± 0.13 mag, one of t he highest among GRBs 
having observed optical afterglows (|Kann et al.ll2010l ). The 
other two models, MW and SMC, are ruled out. Concerning 
the models, a simple power law from NIR to X-rays provides 
a very good result, with /3ox = 1.01 ± 0.05. 

The quality of the photometric set derived for 
GRB 080603A allows us to attempt a full parametric dust 
characterisation using the gene ral parametrisation of t he Lo - 
cal Group extinction laws by iFitzpatrick fc Massal (Il990l) 
(hereafter FM), like with GRB 080607 (jPerlev et al]l201lh . 
In our analysis the parameters ci, C2, and Rv were all tied 
to each other, as for GRB 080607. The other parameters 
that were fixed are 7 to 1 and C4 to 0.6, respectively ac- 
counting for the 2175 A bump width and for the strength of 
the far-UV rise. The free parameters were /? — 0.98 ± 0.04, 
Av,z = 0.57 ± 0.19, Rv = 2.14 ± 0.33, and cg = 1.76 ± 0.66 
(strength of the 2175 A bump) with x^/dof = 19.8/18. 
Indeed, comparing with the corresponding values for the 
LMC2 supersheU, R y = 2.76 ± 0.09 and C3 = 1.46 ± 0.12 
(I Gordon et al ] |2003l l. confirms that the FM model converges 
to values generally consistent with those of the LMC2 pro- 
file. In particular, the 3cr nonzero value of C3 shows that the 
2175 A bump is likely to be present. 

As a further check, we folded the spectral models with 
the filters' transmission curves and iteratively found the new 
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effective frequencies at which the model flux densities were 
the same as the folded ones. We recalculated the best-fit 
model based on the new effective frequencies. This sequence 
was repeated until the best-fitting parameters converged to 
a stable solution. The final effective frequencies changed by 
5% at most from the corresponding nominal values which 
had not been folded with the transmission curves. The best- 
fitting parameters and the possible evidence for the 2175 A 
bump were confirmed with the same confidence. 

Although a broken power-law model cannot be ruled 
out, in this case the break frequency must lie within the 
soft end of the X-ray band. We therefore conclude that the 
break frequency (if any) must lie either outside the optical 
to X-ray range or within the X-ray band itself. 

Figure [7] shows a SED at 4.1 x 10^ s including two radio- 
fiux measurements. Under the assumption that there is no 
break frequency between radio and visible filters apart from 
the peak synchrotron frequency Vm, we tried to model the 
radio points as lying in the Fi, ~ u^^^ p ower-law segment 
for V < in the slow-cooling regime (|Sari et al.l ll998l ') 
with the cooling frequency lying above the X-rays. We ex- 
trapolated the optical flux densities of R and g' , the only 
fllters with observations taken at a comparable epoch, as- 
suming the temporal decay of Section [4. II The X-ray spec- 
trum was also rescaled accordingly. We assumed a negligi- 
ble contribution from the reverse shock. This would allow 
us to constrain at the same time: it turns out to be 
Vrn (1 -I- .z) = (4.4 ± 0.8) X 10^^ Hz in the GRB rest frame 
(dashed-dotted line in Fig. [7|). However, as will be shown in 
Section [5.1.11 this clashes with the observed temporal evo- 
lution of the afterglow, particularly in the radio band: from 
the light curve at 8.46 GHz Vm must have crossed the ra- 
dio band immediately thereafter, around 4.4 x 10'' s. This 
would imply an unreasonably fast decay of u^n- In addition, 
the resulting peak flux of 1.2 mjy is much higher than that 
observed in the radio when crossed it. 

4.3 Optical flash 

Figure [3] clearly shows that the flux densities at optical 
and 7-ray wavelengths during the optical flash that oc- 
curred simultaneously with the last 7-ray pulse are nearly 
equal. The corresponding average spectral index is there- 
fore /3opt-7 ~ 0. Given the considerable amount of dust, 
the dust-corrected optical flux increases by a factor of ~ 7, 
obtained by the LMC2 extinction proflle that best flts the 
SED at the GRB rest-frame frequency of 1.3 x 10^^ Hz, cor- 
responding to the observed R fllter. We point out that the 
dust considered here is that within the host galaxy, since 
the Galactic term had already been removed. Replacing the 
observed optical flux density with the dust-corrected value, 
the average spectral index becomes /3opt-7 = 0.13. The 7- 
ray spectral index during the last pulse is P-y = 0.65 ± 0.2 
(Table [21). 

On the one hand, a simple extrapolation of the 7- 
ray spectrum to optical wavelengths overpredicts the dust- 
corrected optical flux by 2 ± 1 orders of magnitude; this 
has already been observed for other bursts with optic al de- 
tections during the prompt emission (|Yost et al.ll2007l : note 
that these authors adopted a different convention for the 
sign of /?). On the other hand, given the intermediate value 
of lying between the most common values of and 1.3 



expected for the low-energy and high-energy indices (respec- 
tively) of a typical Band function, Ep could lie either within 
or close to the 20-200 keV band with a low-energy index 
close to 0, as we argued in Section 13.11 In this context, 
the dust-corrected fiopt — i value is fully consistent with the 
low-energy index d istribution obs erved in the prompt spec- 
tra of most GRBs (|Kaneko et al. II2OO6) and the optical flux 
density would match the extrapolation of the prompt 7-ray 
spectrum. However, a broadband flat spectrum and a cor- 
respondingly flat electron energy distribution is somewhat 
nonstandard in the synchrotron shock model. 

An interesting possibility proposed in the literature in- 
terprets the optical flash as the result of internal shocks with 
lower velocity irregulariti es at larger radii, as sug gested for 
the flash of GRB 990123 (|Meszaros fc Reej|l99d i. 

Temporal analysis of both proflles adds little informa- 
tion: because of the relatively coarse optical sampling of 
FTN, a correlation between the optical and 7-ray fluxes is 
neither conflrmed nor ruled out. This was tested by inte- 
grating the 7-ray counts in the time windows of the opti- 
cal frames and by comparing the relative variations between 
the two bands. Although the three optical points exhibit the 
same behaviour as the 7-rays, both cases (correlation or lack 
thereof) are compatible with the data within uncertainties. 

Although the measured optical flux can be the low- 
energy extrapolation of the prompt 7-ray emission, in Sec- 
tion 15.21 we test the possibility that the 7-rays are upscat- 
tered photons of the optical flash, as suggested for other 
GRB s to overcome the problems of the synchrotron model 
re.g.. lKumar fc McMahonll2008l . iRacusin et ai]|2008l ). 

4.4 Host galaxy 

We compared the observed host-galaxy SED (Table [3]) with 
a set of spectral syn thesis models by using hyperZ (see 
iBolzonella et al.ll2000l for details). The lack of measurements 
redward of the Balmer break at 4000 A(rest frame) limits 
our possibility of constraining one of the key parameters, the 
total stellar mass. We assumed different extinction proflles 
for modelling the dust extinction produced within the galaxy 
itself; although the best-flt parameters do not change signi f- 
icantly, we adopted the Calzetti law jCalzetti et al.ll200Cll l. 
The result is shown in Figure [6] 

The best-fltting synthetic spectrum belongs to a star- 
burst galaxy with an age of 130 Myr, dust extinction Av = 
0.87 mag, absolute magnitude Mb = —20.7. By adopting 
an LMC extinction proflle, similar parameters are obtained: 
90 Myr age. Ay = 0.77 mag, Mb = -20.7. These values 
are typical of other host galaxies of long-duration GRBs 
(|Savaglio et al.ll2009l : IChen et al.|[2009l ). 



5 DISCUSSION 

GRB 080603A exhibits several interesting properties: (i) the 
end of the 7-ray prompt emission marked by the simultane- 
ous detection of an optical flash, which appears to be a dis- 
tinct component from the emerging afterglow following the 
end of the 7-rays; (ii) an overall achromatic afterglow rise, 
peak, and decay, followed by a break around 10^ s; (iii) an 
accurate SED allowing us to precisely measure the remark- 
able dust content along the sightline to the GRB, clearly 
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Table 6. Spectral energy distribution best-fitting parameters. 



Frequency Range Model Ext. profile /3 Ay z Iz '-'h X'^/^of Prob 

(IQiS Hz) ' (lO^i cm-2) (lO^'^ Hz) (%) 



0.3-3 X 10^ 


POW 


LMC2 


1.01 ±0.05 


0.80 ±0.13 


6 7+'-° 

-l.S 




19.0/19 


46 


0.3-3 X 10^ 


POW 


SMC 


0.92 ±0.04 


0.48 ± 0.07 


r- 0+1.8 
0-°-1.6 




39.9/19 


0.34 


0.3-3 X 10^ 


POW 


MW 


1.04 ±0.04 


0.91 ±0.12 






136/19 


< 10-15 


0.3-3 


POW 


LMC2 


0.8 ±0.7 








0.4/4 


98 


0.3-3 


POW 


SMC 


2 0+"-'^ 

^•'-'-0.6 


< 0.33 






12.1/4 


1.7 


0.3-3 


POW 


MW 


9 9 + 0.2 
^•^-0.4 


< 0.28 






10.5/4 


3.3 


200-3 X 10^ 


POW 




1.3 ±0.3 




"■"-4.6 




16.4/13 


23 



0.3-3x10^ BKN LMC2 0.99 ± 0.07 0.80 ±0.13 S.olg;^ 16.6/18 55 

0.3-3x10^ BKN SMC 0.89 ±0.05 0.48 ± 0.07 T-^ttl ^■'^-I'.o 34.9/18 1.0 
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Figure 6. Observed SED of the host galaxy and the best-fitting 
synthetic spectrum using hyperZ. 



favouring an LMC2 profile at variance with that found for 
most GRBs. In the following we examine these aspects in 
the context of the standard afterglow model. 



5.1 Broadband afterglow modelling 

In the context of t he standard a fterglow model 



ijMeszaros fc Ree3 



Meszaro^ 



19971 : ISari et all 1 19981 : for a review, 
20061) . the power-law piecewise spectra 



see, e.g., 

and Ught curves are interpreted as the result of synchrotron 
emission of a population of shock-accelerated electrons of a 
forward shock ploughing into the surrounding medium. The 
electron energy distribution is assumed to be dN/d'y oc 7"'' 
(7 > 7m); the values of p typically derived from GRB 
afterglow m odelling cluste r arou n d 2.0-2.4 with a sc atter 
of 0.3-0.5 (|Starling et al.1 I2OO8I : ICurran et aO |2010| ). in 

general agreeme nt w ith theoretical expectations (e.g., 

lAchterberg et al.ll200ll:ISpitkovskvll2008l ). 

Following Sari et al. 1 19981 ). let Um and Vc be the syn- 
chrotron injection and cooling frequencies, respectively. The 
observed spectral index /3ox ~ 1-0 (Section 14. 2|) with no 
breaks between optical (fo) and X-ray (i/x) frequencies can 
be explained in two alternative cases: Vni < fo,x < i^c or 
max(i/'ni, i^c) < i^o.x- In the latter case, the electron index is 
p — 2 /3ox ~ 2.0 (although formally consistent with fast cool- 



ing, at late times we can safely assume slow cooling), while 
in the former case p = 2 /3ox + 1 = 3.0 (slow cooling). 



• fm < !^o,x < I'c (P = 3). 

The predicted temporal index a depends on the density 
profile: either oism = 3(p — l)/4 = 1.5 (homogeneous or 
ISM) or Qw = (3p - l)/4 = 2.0 (wind). Both decay val- 
ues are significantly steeper than the observed Q2 ~ 1.0, 
respectively, by Aqism = 0.5 and Aa^ = 1-0. Even as- 
suming a more general density profile, n(r) « r"", the 
expected decay is steeper than ISM for every value of s. 
This was also the case for a number of other bursts whose 
spectral and temporal indices could not fulfil any closure 
relation (e.g.. iMelandri et al. 2008). A possible way to ex- 
plain a shal lower decay is energy i njection refreshing the 
blast wave (iRees fc Mfezaroslll99^: ISari fc Meszarod I2OO0I : 
iGranot et al.ll2003l : IMelandri et al.ll2009l ). as was also pro- 
posed to explain the shallow-decay phase in early X-ray af- 
terglows (Zhang et al. 200J3). 

Let E{t) be the fireball energy as a function of the ob- 
served time t, so that E{t) oc t"^. Assuming negligible ra- 
diative losses, the expected decay index change with respect 
to no injection, A aei, is e(p + 3)/4 (ISM) and e(p + l)/4 
(wind) for u < Uc (|Panaitescul |2005| ) . The observed values 
imply e = 1/3 (ISM) and e = 1 (wind). The energy bud- 
get may be problematic, as in the afterglow phase from 10^ 
to 10^ s the required injected energy would be a factor of 
10 0^/^ 5 (100) larg er for ISM (wind). Using the notation 
of IZhang et al.l (I2OO6 ). these are equivalent to 5 = 2/3 and 
q = 0, respectively, for an injection luminosity L{t) oc 

At the end of the injection, the power-law decay is ex- 
pected to steepen by Aoei so that the decay should resume 
the no-injection values, aisM or a„. Indeed, the final decay 
Q3 — 1.71q'3 compatible with both values. In this case, 
the final break around 10^ s would mark the end of the 
energy-injection process. 

• max(i^m, i^c) < J^o,x (p = 2). 

Above the cooling frequency the emission does not depend 
on the ambient density, so for any s (0 < s < 3), a = 
(3p — 2)/4 = 1.0. This fully agrees with the observed value 
of 02 and there is no need to invoke any additional processes 
such as energy injection. Another asset of this possibility is 
that the final power-law decay index, 03 ~ 2, nicely supports 
the jet-break interpretation, being a = p. 
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5.1.1 Radio afterglow 

The peak frequency of the synchrotron spectrum crossed 
the radio bands between 5 and 12 days: this is suggested by 
the peak in the light curve at fradio,2 = 8.46 GHz and by 
the change in the spectral slope between !-'radio,i = 4.86 GHz 
and z^radio,2. In particular, the slope of the radio spectrum 
changes from positive to negative: the first 2-channel ra- 
dio SED at 5 days is fit with /?radio = — 1.3^2' J; while the 
second SED around 12 days gives /3radio = l-Sl^ g. At low 
frequencies the flux is expe cted to rise as t ^^'^ and to decay 
as f-3(p-i)/4 for an ISM (|Sari et al.llT998l ). We fitted the 
observed 8.46 GHz radio curve under these assumptions in 
either case considered above (i.e., either p = 3 or p = 2) 
and found the radio peak time iradio.p = 4.4lQ y x 10^ s 
and the peak normalisation -Fradio.p = 200^40° from 
the fit. This could be considerably different, when other 
temporal behaviours are considered; in particular, a steeper 
rise followed by a steeper decay could be compatible with 
-Fradio.p > 200 /iJy. This possibility is not naturally ex- 
plained within the standard afterglow model, unless the //m 
passage and the jet break happened almost at the same time 
by chance, implying a decay. 



5.1.2 Afterglow onset 

The rise experienced by the afterglow in the BVRi filters 
rules out the passage of a typical synchrotron frequency 
through the observed wavelengths as the possible cause be- 
cause of both the steepness and the lack of chromatic evolu- 
tion of the rise itself. The possibility of an afterglow emerg- 
ing from a wind surrounding a massive progenitor, with the 
optical rise being due to the progressively decreasing dust 
extinction (|Rvkoff et al.|[20o3 ). is excluded by the lack of 
chromatic evolution , as observed in man y other cases (e.g., 
iGuidorzi et ai]|2009l : iKriihler et a"l]|2009l ). 

Here we consider the possibility, disc ussed in several 
analogous cases (e.g., Molinari et"alll2007l '). that the broad 
optical peak tp = ISTSljso ^ niarks the deceleration of the 
fireball and the afterglow onset. The duration of the 7-ray 
burst itself (~ 200 s) is much sh orter than the peak time, 
as expected in the thin-shell case (|Sar il ll997f ). The observed 
steep rise, qi = —3.6 ± 1.7, rules out the wind environment 
for which a shallower rise (~ —0.5) is required. Depending on 
whether it is !^ < fc or > Vc, a rise index of —3 or —2 is ex- 



pected for an ISM (|jin fc Fanll2007l :[ Panaitescu fc Vestrandl 
l2008l 'l. both compatible with observations. 

In this context, from the afterglow peak time we can 
estimate the initial bulk Lorentz factor To as being ap- 
proximately twice as large as its value at the peak time 
jSari fc Piranlll999l : [Moiinari et al.ll2007^ : 



To 



2r(tp) = 2 



3giso(l-Hz)^ 
32TTnmpC^rjtp 



1/8 



(130±20)n~^''**(3) 



We assumed standard values for the energy-conversion ef- 
ficiency, rj^ = 0.2 and for the particle density of the cir- 
cumburst environment n = ng cm~^. Equation ^ holds for 
an ISM environment; we do not consider the wind case, be- 
cause of the incompatible steepness of the rise. Such a bulk 
Lorentz fac tor lies within the distribution found f or other 
GRBs fe.g.. lLiang et al.ll2010l:lMelandri et"ai]|2010l '). 

Following IZou fc PiranI l|2010l v" we set an upper limit to 



To from the prompt 7-ray light curve, thanks to the pres- 
ence of a quiescent time between the two pulses. The idea 
behind this is that while the prompt 7-rays are being pro- 
duced through internal shocks, the outermost shell begins 
sweeping up the surrounding medium; as a consequence, 
a forward shock should appear in soft 7-rays as a smooth 
and continuously increasing emission. This constraint is only 
suitable for bursts with a relatively short pulse followed by 
either a quiescent time or a deep trough, because otherwise 
the external shock could have significantly decelerated dur- 
ing the first pulse. From Figure |4l at 90 s we estimate a 3a 
upper limit of 3 fiJy at the mean energy of 84 keV, which 
corresponds to 0.3 x 10~^* erg cm~^ s~^ Hz~^. From their 
Equation (5) we derive the following upper limits to Fq: 

Fo < 150n-i/«.-/f/,.-Zf (l + y)^/^« (p = 3) (4) 
Fo < 220n-^/^:L'',%(l + r)^''« (P = 2.01), (5) 



where we used p — 2.01 instead of p = 2 as it is applica- 
ble only for p > 2. A value p < 2 would imply electron 
energy divergence, so a proper formulation is required; this 
is beyond the scope of the paper, es = €3,-1 x 10~^ and 
Ee = ee,-i/2 X 10"^''^ are the equipartition factors for the 
magnetic and the electron energy densities, respectively, and 
Y is the Compton parameter for synchrotron self-Compton 
scattering. These upper limits are remarkably consistent 
with the estimate derived in Equation ([3}, especially be- 
cause a number of GRBs were found to exceed these values. 

In the context of peaks interpr eted in terms of outflow 
decele ration and afterglow onset, IPanaitescu fc Vestrandl 
(|201ll ) proposed that fast-rising optical afterglows are likely 
caused by an impulsive ejecta release with a narrow distri- 
bution of Lorentz factors after the GRB itself, in contrast 
to an extended release or a broad range of Lorentz factors 
more suitable to explain the slow-rising/plateau afterglows. 
The motivation of this interpretation resides in the differ- 
ent correlations between peak time and peak flux found for 
each class: Fp oc tp"^ and Fp cc tp^^'^ for fast-rising and 
plateau afterglows, respectively. Following their guideline, 
when we move GRB 080603 A to a common redshift of 2 = 2, 
its dust-corrected 7?-band peak flux is 0.8 mjy at a peak 
time of 1750 s (at z = 2) — that is, it lies in the region of 
the Fp-tp pla ne where the two correlations cross each other 
(see Fig. 1 of 'Panaitescu fc Vestr andll201ll ). Because of the 
~ rise, GRB 080603 A belongs to the peaky afterglow 
class. Interestingly, _Eiso and Fp of peaky GRBs are shown 
to correlate more tightly than those of plateau GRBs; indeed 
GRB 080603A lies very clo se to the best-fitting power-law 
relation shown in Fig. 2 of IPanaitescu fc Vestrandl l|201ll ). 
in agreement with its being a member of the peaky class. 
The tighter connection between the 7-ray released energy 
and the afterglow peak flux for the fast-rising GRBs may 
support the impulsive ejecta release interpretation. 

In the afterglow the presence of a single peak fol- 
lowed by a ~ t~^ decay qualifies GRB 080603 A as a Type 
HI me mber according t o the classification by IZhang et al.l 
(1 200^ and |jin fc Fan! (120071 ). At variance with pre- Swift 
expectations mainly based on the case of GRB 990123 
(jAkerlof et al.l [T999I '). most GRBs show no evidence for a 
short-lived reve rse-shoc k peak at early times (|M"undell et al.l 
| 2007l:IOates et a l. 2009; Rvkoff et al.ll2009l : lKann et al.ll2010l : 
Melandri et al.|[2010 ) and GRB 080603A is no exception in 
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this respect: the prompt optical flash cannot be reverse- 
shock emission because of the gap between the flash and 
the afterglow onset. A way to circumvent this problem is 
that the outer layer of the outflow has much higher Lorentz 
factor compared to the bulk part of the flow; the outer layer 
might produce an optical reverse-shock emission well before 
the onset of the afterglow, which is determined by the main 
part of the flow. However, the rather narrow spike seems 
to disfavour this model. Among the possible explanations, 
such as either a high or a low magn etic energy density in 
the ejecta (e.g.. ICombo c et al.^l 2008l). here we consider the 
low-frequency model "([NIelandri et al.l I2OI0I ) : at the shock- 
crossing time, marked by the peak, both injection frequen- 
cies of forward and reverse shock s, i/m.f and Vm r (res pec- 
tively), lie below the optical band ijMundell et al.il2007l ). 

Assuming the same microphysical parameters in both 
shocks, the relation between the spectral characteristics of 



(6) 



where Fmax is the peak flux in the frequency domain at a 
given time, in this case at tp, different from the peak flux in 
the time domain at a given frequency, denoted with F-p. We 
discuss the implications in the two cases considered above. 

• I^m < Vo,^ < i^c {p = 3). 

This requires i/ni,f ^ 1^0 and i^c > J'x at t = tp, where 
Uo = 5 X 10^* Hz and ~ 10^* Hz. As discussed above, 
the observed temporal decay requires energy injection with 
g = 2/3 up to the late-time break. Temporal evolution of 
the characteristic frequencies and peak of the forward shock 



(9+2)/2 



t 



-4/3 



OC t 



(q-2)/2 



t 



-2/3 



and 



is Vm,! CC t 

i^max.f OC t^-' = t^^^ (|Zhang et al.ll2006l ). Interpreting the 
radio flux at tradio = 4.1 x 10^ s as mostly due to the forward 
shock, using the result of the corresponding SED fitting dis- 
cussed in Section r4.2l and shown in Figure [7] (dashed-dotted 
line), we can estimate the time at which i/^^f crossed the op- 
tical bands: tradio (i^o/j^m.f (tradio))"^^* (tradioAb2)^'^* = 7 X 

10^ s. This is derived from the temporal dependence of 
i^ni.f OC t~'^^^ as long as energy injection goes on (so for 
t < th2), and fm.f OC t~^^^ (for t > th2)- The absence of any 
break or spectral evolution in the multi-filter light curve of 
FigureUrules out any such passage at this time. Even worse, 
as discussed in Section 14.21 the derived value of !/m,f from 
Figure [7] is incompatible with its passage through the radio 
bands at tradio.p as observed from the radio curve (Section 
I5.1.1|) . because it implies a too rapid decay of z/m.f. 

The expected contribution to the radio flux from the re- 
verse shock exceeds the observed values by a factor of 10^: 
at the peak time, Fmax.r ~ 80 mjy. This is found from 



Equation ^ and from 



F, 



p,uiicxt 



= 1200 /iJy, 



where we used the observed flux density in the R band 
at the peak time, Fp = 180 ± 6 fJ.3y; removing the in- 
trinsic dust extinction of Av,z ~ 0.8 mag, the unextin- 
guished value is -Fp.unext = 1200 /iJy. From the decelera- 
tion time to the end of the energy injection, i t is -Fmax.r ~ 
Fmax .f r OC t^-i-(2+i)/» = ^0 (Table [3 and Zhang et al.' 
200(: ). From that time to ^radioi it is Fmax.r oc t^^ (Kobavashi 



2000l : IZhang et al] 120031 '). It is therefore f;nax,r(tradio) 
80 (tradio/ib2)~^ ~ 20 mJy. From Equation ((B)) z/m.r 



1.2 X 10 Hz at the peak time. Using the tem- 




jglO ^q12 ^gl3 ^q14 ^q15 ^q16 ^q17 ^q18 

Vrest [Hz] 

Figure 7. Rest-frame spectral energy distribution at 4.1 X 10^ s. 
Two alternative models are shown, depending on whether lies 
below the optical (p = 2; solid line) or above X-rays (p = 3; 
dashed-dotted line). Characteristic frequencies are reported. Both 
models are extinguished with an LMC2 dust profile. The dashed 
line is the same as the solid line, removed of dust extinction. 



I^m,r(iradio) ~ t'm.r (ib2/ip) '^^^ (tradio/ib2) '^^'^ ~ 1 GHz, 

where we used v^n,! oc t~^^^ for t < tb2 (Table [7| and 
Um,r OC for t > tb2 (|Kobavashill2000l ) . This means that 

around the time of observations, Um,! should have crossed 
the radio band. The expected flux should therefore be com- 
parable with fmax,r(tradio) ~ 20 mJy. Although the passage 
of the peak synchrotron (of the reverse shock) through the 
radio band is indeed compatible with the radio light curve, 
the observed flux is about two orders of magnitude smaller 
than expected. Although in principle self-absorption could 
explain this, in the following we show that here it appears 
unlikely. Another possibility to suppress the reverse-shock 
radio emission at late times is assuming a hot flow (i.e., 
Poynting-flux dominated) . 

Equation (|6]) is valid as long as we consider energy injec- 
tion to the shock by a cold flow (i.e., injection of kinetic en- 
ergy). More generally, assuming the same velocity and pres- 
sure in the FS and RS regions, the ratio Fmax,r/-Fmax,f is pro- 
portional to the corresponding ratio of the number of elec- 
trons in the two shock regions. Considering an extreme case 
of no new electron injection in the reverse-shock region, it is 

F^ax.r/F^ax.f « A^e,r/A^e,f ~ r"^ ^ t-3(2-9)/4 ^ ^-1^ 



3(2-<j)/2 



r / -fma; 

larly, !/m,r/fm,f « (7e,r/7o,f)^ ~ {Pi / Pif ~ r" t 

t^; i/c,r ~ I'cf ~ These might be more un- 

certain compared to the cold-flow results, where we have as- 
sumed that the FS and RS regions have comparable widths. 
• max(i^m, J'c) < J^o,x (p = 2). 



This requires u^-^s < Vo and Uc < 



at t = tp. No en- 



poral scaling of !^m,r, its value at fradio is found to be 



ergy injection is required to explain the spectral and tem- 
poral properties of the afterglow. The characteristic syn- 
chrotron frequencies can b e expressed by (jSari et al.., 19981 : 
iKobavashi fc ZhandlioOSal ) 

/, \ no inl2 -3/2^1-1/2 1/2 -l*-l/2,-l/2 „ /r,\ 
I/c(tp) ~ 6.8 X 10 -E^52 %.2% C ip,3 Hz, (7) 

l-t \ ^ T '7 inlS 1/2 2 Tpl/2 -1/2^1/2,-3/2 „ 

Vi-n,i(tp) ^ 1.7 X 10 e^' e^E^ij r^g^ C tpj Hz, (8) 

where rjo.2 = r?/0.2, tp,3 = tp/10^ s, Eb2 = -Eiso/10^^ erg, 
and = (1 -f 2;)/2.687. In particular, it is ~ 3.6 x 
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Shock L oc E{> 7) oc 7-^+1 



FS -(g + 2)/2 (g-2)/2 (1 - g) -12/(7 + s) -2(s + l)/(7 + s) 3(s - l)/2(7 + s) 

RS -(g + 2)/4 (i2-2)/2 3(2 - 3g)/8 -6/(7 + s) -2(s + l)/(7 + s) 3(s - 2)/2(7 + s) 



Table 7. Temporal exponent of injection and cooling frequencies, as well as of the maximum flux density for the forward shock (FS) 
and the reverse shock (RS) in the case of continuous energy injection through an ISM. The two formalisms, the luminosity as a function 
of time on the left- hand side, and the energy distribution o f the shells as a function of the bulk Lorentz factor on the right-hand 
side, are equivalent jSari fc Meszarosll2000l : IZhang et alll2006h . The two parameters q and s are related by s = (10 — 7q) /(q + 2) and 
q = (10 — 2s)/(7 + s). The impulsive case (i.e., no continuous injection) corresponds to s = q = 1. 
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Figure 8. Top panel: rest-frame SED at 4.1 X 10^ s. The solid line 
shows the synchrotron spectrum with an LMC2 dust-extinction 
profile for the p = 2 case. The corresponding injection and cooling 
frequencies are indicated. The dashed line shows the same unex- 
tinguished model. Bottom panel: rest-frame SED at 1.1 X 10^ s. 
The data are consistent with being taken after the passage of 
through the radio band. 



-3/2 



- 1/2 



Hz, 
2 Hz. 



1.3 X 10^* ei/^ee Hz, and 



The requirements on the characteristic frequencies at the 



peak time become ej/^te 4 x 10~' 
A possible solution is es « 10"^, 
which gives Vc ~ 3.6 x 10^'* Hz and t'm.f 
The much lower i/c than the previous case would be mainly 
due to a denser environment. 

The expected peak flux of the forward shock 
at the peak frequency j/m.f is then -Fmax,f(ip) ~ 
(j^m,f/z^c)-'''-''/'(/^c/j/o)-''/'-Fp,unext ~ 3 mjy. The 
contribution of the reverse shock to the optical peak 
luminosity in the R band is comparable: from Equation ((6| 



and e~^^^n~' < 140. 
f» 3 X 10"^ no ~ 10, 
1.2 X 10^'' Hz. 



the reverse shock peaks at Vm.i ~ I'm.f/r « 2.4 x 10 Hz 
with -fmax.r ~ T Fjnax.f which scales at fm.f by the 
factor (/^m.f/i^m.r)"*''"^''^^ = T^V The net result is 

At t = tradio, the expected luminosity is dominated by the 
forward shock. From the frequency and flux-density scalings 

I'm,! oc t"^/^, Vc,f oc t~^'^^ Fmax.f OC t° Vm.r OC t'^^'^, Vc.r OC 

tZHI' and Fmax,r OC t'^ (jzhang et al.ll2003l : iMundell et all 

[2003), we have !-'m,f (tradio) « 3 X 10^° Hz, t'm,r(tradio) « 7 X 

10^ Hz. This means that the peak of the reverse shock has al- 
ready crossed the radio band, while that of the forward shock 
has not yet. The two expected flux densities at 8.46 GHz 
are F/(iradio) ~ 3 (i^radio/j'm,f (tradio))^''^ « 2 mJy and 

-Fr(iradio) ~ 3 (l^radio/l^m.r (iradio) ) ~ ^'''^ (tradio/tp) ~ ^ mJy 

« 0.3 /iJy, respectively. Thus, given that the expected flux 
from the reverse shock is 3 orders of magnitude smaller than 
observed, while that of the forward shock differs by a factor 
of few, hereafter we focus on the latter. 

The solid line of Figure [7] shows the result of fitting the 
broadband SED at tradio with both i/m.f and fcj below the 
optical and having fixed p = 2 and the rest-frame dust ex- 
tinction. The free parameters are the break frequencies as 
well as the normalisation. We found Um,! (1 + ^) — l-T^g^ x 
10" Hz and i/c,f (1 + z) = Stl x lO" Hz (xVdof = 20.9/17). 
Not only does the value for agree with expectations, 
but also the peak fiux of 400 fiJy resulting from the flt 
is nearly compatible with the value derived from the radio 
light curve (Section lS.l.ip . The downside is the value for i^c.f: 
being so close to the optical band, it must have crossed it at 
some time t < tradio- This would imply /3o = (p— l)/2 = 0.5 
for t < tradio, which is clearly not true. A possible solution 
could be to assume a cooling frequency that increases with 
time, but given that we ruled out a wind environment, this 
option is not acceptable either. Another problem concerns 
the crossing time of i^c.f through the radio band, observed 
immediately afterward at iradio.p ~ l.liradio (Section lS.l.lf) . 

We tried to decrease the latter problem by flxing v^j to 
the value expected from the time it crossed the radio band 

that is, by imposing i/m.f = J^radio (tradioAradio,p)~^''^ = 

9.4 GHz. We allowed all the remaining parameters to vary. 
The result is shown in the top panel of Figure [H] The best- 
fltting parameters are Ucj {1 + z) — (1.3 ± 0.1) x 10^^ Hz, 
/3ox = 0.92 ±0.04 and Av,z = 0.56 ±0.10 mag. Although the 
spectral index is slightly harder than that found at previous 
epochs, both the normalisation and the slope of the radio- 
to-optical spectrum flt in the expected broadband modelling 
well. This explains both the spectrum and light curve of the 
radio observations. The bottom panel of Figure |5] shows the 
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SED at t = 1.1 X 10® s, when v-^^f has just crossed the radio 
band. Apart from the same issue with i^c.f already men- 
tioned, this is in remarkable agreement with expectations: 
y^ f (1 + 2) = (1.3 ± 0.1) X 10^5 Hz, /3ox = 0.98 ± 0.05, and 
Av,z = 0.72±0.10 mag. Yet, the inferred temporal evolution 
of the cooling frequency remains an issue. In this case, this 
would lie within the optical bands and the aforementioned 
argument still applies. 

Overall, the p — 2 case works better than p — 3 and can 
account for more observed properties of the broadband af- 
terglow evolution. Still, the derived evolution of the cooling 
frequency with time conflicts with a homogeneous environ- 
ment, the only one compatible with the data. 



5.1.3 Self-absorption 

So far we assumed a negligible effect due to self-absorption 
in the observed radio flux. Should the radio flux be self- 
absorbed, from Figure[7|both cases could be compatible with 
having a high value of Fmax. In particular, for the p — 2 case 
the temporal evolution of fc,f would no longer be an issue, 
because it could lie well below Vo Si.t t — tradio. To signifi- 
cantly suppress the radio flux, one should require Vc <S i^o 
already at the peak time; however, this would imply a much 
larger optical luminosity and an unusually high energy bud- 
get. This seems very unlikely, given that the optical lumi- 
nosity of GRB 080603A, corrected for the dust extinction, 
already lies in the mid-to-bright end of the observed optical 
afterglow distribution (Section 15.31 and Figure [TT]). 

Although at lower frequencies and early times self- 
absorption can significantly suppress the flux, this does 
not explain these observations unless one makes extreme 
assumptions. A simple estimate of the maximum flux is 
that of a black body with the forward-sho c k temperature 
dSari fc Piranlll999l : iKobavashi fc Sarill2000l : iMundell et al.l 
I2OO7I ). which at the peak time in the optical {t = tp) is given 

by 



F„.BBitp) fa 71(1 + z)i''^eempr'^ 



500 



^radio 



R -2 

L 

1/2 
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(10) 



where R± ~ A.QTctp is the observed fireball size and the 
dependence on no is inherited from Equation Q . The value 
derived from Equation (|10|) initially increases as ~ t^^^, and 
steepens to ^ f^^* after z^ni,f crosses the observed frequency; 
thus, at the time of radio observations the black-body flux- 
density limit expressed by Equation (|10|l increases by a fac- 
tor of ~ 20. This can hardly explain the observed radio flux 
unless one assumes ~ 10~* and/or a high-density envi- 
ronment (n « 400 cm~^). This value for ee would imply 
that Um is much below the optical bands (t^m.f ~ 10*~^ Hz) 
at the peak time and, consequently, an unreasonably high 
value for Fmax.f. 



5.1.4 Off-axis jet 

In the off-axis jet interpretation, the afterglow rise and peak 
do not mark the fireball deceleration, but are the result of a 
geometric effect due to an observer angle, Sobs, being larger 
than the jet opening angle, 9j: as the jet decelerates, the 



relativistic beaming angle l/F progressively increases, re- 
sulting in a rising light curve as seen from the off-axis ob- 
server, provided that the jet structure is such that the en- 
ergy and F drop sharply at 9 > 9j. The peak is reached 
when F (Sobs — Sj) ~ 1, as the sightline enters the beam- 
ing cone at the edge of the jet; finally, the decay asymp- 
totically approaches the light curve for an on-axis observer 
(jGranot et al.l|2005l ). 

We can estimate Sobs from the peak time as follows: 
Sobs — Sj ~ 1/F(tp). Interpreting the late-time break as be- 
ing due to the jet, it is Sj « 1 /F(ij ) , from which we estimated 
Sj ~ 5?7 (Section l5.3p . Assuming that the deceleration of the 
jet occurs before the peak, we can assume the temporal scal- 
ing of t he bulk Lorentz factor a s F cx t~^^^ for an adiabatic 
cooling l|Rees fc Meszaroj[l99^ : 



r(tj 

F(tp 
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0.2, 



(11) 



so (Sobs-Sj) « 1°. 

Given the steep rise observed for GRB 080603 A, we 
conveniently consider the case of a sharp-edged, homoge- 
neous jet seen at Sobs > Sj. The observed peak energy 
of the prompt emission, Ep, falls off as fe^, where b — 
F ( Sobs — Sj), while the observed energy falls off rapidly as 
&® (jGranot et al.ll2002h . From Equation dTT]), b ^ 2Fioo, 
implying the following on-axis values: £'p(Sobs = 0) ~ 
eOOFfoo keV and, and £;iso(Sobs = 0) lO''* Ffoo erg; here, 
F is the initial value of the bulk Lorentz factor. These esti- 
mates are only illustrative and show that, in principle, this 
scenario could work for reasonable values of the physical 
parameters in the case of a uniform sharp-edged jet. 

However, more realistically the jet cannot be exactly 
uniform with very sharp edges; in particular, F is expected 
to be lower at the edge and higher at the jet core. Hy- 
drodynamic simulations have shown that, particularly for 
1 ^ Sobs/Sj < 2, very shallow rises or decays of the early- 
afterglow light curves are expected for realistic je t structure 



afterglow light curves are expected tor realistic ,ie t structur 
and dynamics (|Granot et "aril2002l ; Isichler fc Gra not 200t| 
Alternatively, structured jets l|Zhang fc Meszarosll200' 
with most of the energy concentrat ed in the core can repro- 
duce steep rises in the light curves (IKumar fc Granotll2003l : 
iGranot et al.ll2005l: lEichler fc Granotll2006D. 

Following Panaitescu fc Vestrandl l|2008l ). assuming an 



angular profile for energy as £{6) oc (S/Sj)~', & ^ t rise 
would be obtained for q = 4 and Sobs/Sj ~ 1.5 ^ 2.5. 
GRB 080603A lies within the la region of the correlation be- 
tween peak time and peak fiux for fast-rising afterglows, also 
discussed by Panaitescu fc Vcstraiid (2011) f Section 15.1.21 
see Fig. 2 of fPanaitescu fc Vestrandll2008l ). 

The scenario of a viewing angle slightly larger than 
the jet angle seems to be consistent with the nature of a 
typical GRB, in contrast to the classes of spectrally softer 
events, such as the so-call ed X-ray rich GRBs (XRRs) or 
X-ray hashes (XRFs; e.g., ISakamoto et al. I l2008l and refer- 
ences therein) . Larger viewing angles are expected to be as- 
sociated with softer observed events ijYamazaki et al.l I2OO2I : 
iGranot et al]|2002l : iGranot et al.ll2005^ . and the off-axis in- 
terpr etation for events of this type appears to be favoured 
fe.g.. iGuidorzi et ahllioogl) . In this context, GRB 080603A 
could represent a soft /intermediate classical GRB with a 
typical jet opening angle and viewed with a comparable 
viewing angle. 
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5.2 Prompt optical/7-ray emission: inverse 
Compton? 

Before the onset of the afterglow marked by the steep rise 
at the end of the 7-ray prompt emission, the optical flux de- 
tected simultaneously with the second and last 7-ray pulse 
is unlikely to be synchrotron radiation of the shocked ISM. 
As noted in Section 14.31 the spectral index measured dur- 
ing the prompt emission in the 7-ray band is likely to be an 
intermediate value between the typical low-energy and high- 
energy photon indices of a Band function. In this respect, 
the optical emission could be consistent with the extrapo- 
lation of the 7-ray spectrum down to the optical band. A 
cross-correlation study between the optical and 7-ray pro- 
files would certainly settle this issue; however, in practice 
this is not possible because of the coarse optical coverage, 
which gave only three points separated by gaps in between 
(Fig. ID). 

The variety of observed behaviours in other GRBs is 
rich: the prompt optical was observed to be uncorrelated 
with the ong oi ng hi gh-energy emission for GRB 990123 
Akerlof et all lll999f ) (e.g., see also GRB 060111B, 



Klotz et al 



Stratta et al.1 l2009l : GRB 080607, 

Perlev et al.l I2OIII ). whereas a strong corr elation was ob- 
served, for example, for GRB 050820A (|Vestrand et al.l 



2006 ) , superposed on the onset of the afterglow ( Cenko et al.l 
20061 ). Similar cases of some degree of correlation between 



the 7-ray prompt and o ptical emiss i ons ar e GRB 041219A 
' Vestrand etahl jgOOSi : [Blake et al.1 120051 ') . GRB 060526 



iThoiie et al.l I2OI0I'). an d GRB 080319B jRacusin et al.1 
I2OO8I : iBeskin et al.l l2oiol l. One of the most common cases 
is that the prompt optical observations, typically starting 
during the final part of the 7-ray emission, suggest the tran- 
sition from the inner-engine activity to the multi-band after - 
glow onset; se e, for instance, GR B 051111 dYost et al ]l2007h 
GRB 081008 (lYuan et al.lboiol'). GRB 081 126 (iKlotz et al.1 
l2009h . and GRB 080928 (|Rossi et al.ll201lh . In one case, a 
strong optical flare incompatible with an exte rnal-shock ori- 
gin w Eis observed before the afterglow onset (|Greiner et al.l 
l2009h . In some other cases, the optical profile is dominated 
by the onset of the extern al shock of the e jecta through 
the ISM (e.g ., GRB 080810, 'Page et al.l l2"009l : GRB 061007, 
iRvkoffet all [2009: MundcU et al. 20ol). 

Evidence has also been reported for a sizable temporal 
lag of a few seconds betw e en the optical and high-energy 
profil es (jKlotz et al.l |2009| : iRossi et al.1 12OIII: iBeskin et al] 
I2OI0I ). This potentially represents a strong clue to explain 
the prompt-emission mechanism. 

Comparison of the optical and 7-ray fluxes is also useful 
for establishing the possible link: while some GRBs have an 
optical-to-7 spectral index /3opt-7 compatible with /3-y, as 
can be the case for GRB 080603 A, other events show an 
excess of optical emission with respect to the extrapolation 
of the high-energy spectr um: up to 10^ time s larger, as in 
the case of GRB 080319B (jRacusin et al.ll2008l '). By contrast, 
there are also cases in which the optical flux lies below the 
high-energy spectrum extrapolation, such as GRB 050401 
ijRvkoff et al.1 12OO5I : [Vost et al.1 12007| ). The latter situation 
does not necessarily imply different origins or mechanisms 
for optical and high-energy emissions, but could merely be 
due to either dust extinction (as is at least partially the 



case for GRB 050401, iKann et al.ll2010l ) or a synchrotron 
spectrum peaking between the two energy ranges. 

GRB 080603A is an example where both components 
(internal activity and afterglow onset) are clearly temporally 
separated. Although no firm conclusion can be drawn on the 
possible existence of a temporal lag between op tical and 7- 
ray p hotons, this GRB resembles GRB 081126 (jKlotz et all 
'2009): both 7-ray proflles consist of two disjoint FRED-like 
pulses, the last of which is observed simultaneously with an 
optical flash, f ollowed by the afte rglow onset. 

Following [Piran et al.l (|2009l ). we tested whether the 7- 
ray prompt emission of GRB 080603A can be explained in 
terms of an IC process by a population of relativistic elec- 
trons on low-energy seed photons. The same electrons would 
also upscatter the 7-ray photons to GeV-TeV energies. To 
avoid the energy crisis (i.e., when most of the energy is re- 
leased in the GeV-TeV range owing to a large Y parameter) , 
we used the simultaneous optical and 7-ray flux densities to 
constrain the bulk Lorentz factor F. Let 7e be the Lorentz 
factor of electrons within the fluid rest frame, i^l the (GRB 
rest-frame) peak frequency of the lower spectral component 
(i.e., of the seed vF^ spectrum), and Fl the corresponding 
peak flux. Also, let Uopt = Uo{l + z) = 1.34 x 10^^ Hz be the 
rest-frame frequency corresponding to the observed optical 
band. Fl can then be expressed as 



(12) 



where Fopt = 160 /iJy is the unextinguished flux of the op- 
tical flash. Two possible cases are considered: for the UV 
(IR) solution, corresponding to vi^ > Vopt {v-l < I'opt), we 
assume three possible values for the spectral index: 13 = 0, 
—0.5, and —1 (/3 = 1, 1.5, and 2). The Compton parameter 
Yl in the flrst IC scattering is 

-(1-/3) . / \ -(1-/3) 
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where hv^ = 84 (H~ 2) = 226 keV and = 41 ^Jy (Ta- 
ble [l|. The first-order IC scattering is not in the Klein- 
Nishina (KN) regime, so it is v-y/vj^ = 7^, or, equivalently, 



7. = i^Y' r-^Y' 

\Vopt/ \ VI, J 
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The second-order IC scattering might be in the KN regime, 
so Yh is 
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where ^ = (7c/r) hv-y/rrioC? is the correction factor in the 
KN regime. The energy of the upscattered photons is 



,.«.,CeV(i)%.,..(,IlS^). 

and the corresponding Compton parameter Fh is 
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The energy-crisis problem mainly resides in the large 
value for Yh. Figure [9] shows Yh as a function of vi, for the 
values of the spectral index considered above and assuming 
for r the value we derived in Equation There are two 
solutions for which the released energy is not an issue: the 
UV solution requires /3 ^ — 1 and vi, > 30 Vopt , equivalent 
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VL/Vopt 

Figure 9. Yh as a function of vi^/vopt for F = 130 and (from top 
to bottom) /9 = 0, —0.5, and —1 for i^l > Vopt, and (from top to 
botto m) 13 = 1, 1.5, 2 for vj^ < i^opt (adapted from IPiran et aP 
l2009t) . 



to 7c < 37 (Eq. H14[)). This solution is characterised by a 
neghgible KN suppression. The total energy is at least ?>E~,, 
given that Fh ~ ^l- Other problems with the UV solutions 
are (i) a low efficiency due to the low 7c, because protons 
would carry at least a factor of mp/7e me more energy than 
electrons, unless one requires pair loading; and (ii) at such 
low values of 7e, increases and its flux _Fl is limited not 
only by -Fopt, but also by the prompt soft X-ray observa- 
tions .^WliBneve^avail^ these data rule out the UV solu- 
tion (|Piran et al.ll2009l ). Although for GRB 080603A there 
are no prompt X-ray observations, the UV solution appears 
contrived. 

The IR solution requires /3 ^ 2 and vi^ < O.Oli/opt, 
equivalent to 7c > 2000. This solution is characterised by a 
high value for 11, and a relatively small value for Yh, because 
of the KN suppression. This becomes important for i/l < 
QAGvo-pt, as shown by the break in Figure [9] On the other 
hand, for very low i^l /I'opt the expected _Fl also increases for 
a given _Fopt and self-absorption can represent an issue. To 
avoid this, one has to require -F'sa(i^L) > -Fl, where Fsb.{v-l) is 
the black-body flux for a local temperature kT « r7e ma(? , 
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where R is the radius of the source and dL = 12.74 Gpc is the 
luminosity distance of GRB 080603A. From this requirement 
on _Fl we can constrain F for vi^ < Vo-pt- 



F < 



1 



For 
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We conservatively assume R = 10^^ cm. The corresponding 
allowed F-7e phase space is shown in Figure [10] in the IR- 
solution domain, enclosed by the decreasing curves on the 
right-hand side. Combining this with the Yh ^ 1 regions 
and with the corresponding values for the spectral index 
shows that the value measured for F in Equation ^ does 
not overlap with any total allowed region (darkest areas). 

We conclude that the prompt optical and 7-ray data are 
not compatible with an IC origin for the latter as a result of 
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Figure 10. The allowed (shaded) phase space characterised by 
^ 1 (from bottom to top) for /3 = 0, —0.5, and —1 for i'l > 
i^opt and (from bottom to top) /3 = 1, 1.5, 2 for i/L < fopt- The two 
decreasing functions in the IR solution mark the self-absorption 
limits for /3 = 2 (grey area; lower 7e range) and /3 = 1 (higher 7e 
range), respectively. The two darkest areas partially overlapping 
each other are the intersection of the allowed regions for f} = 2 
and for /3 = 1, respectively. The interval F = 130 ± 20 (estimated 
from the peak in the afterglow light curve) is highlighted. 



upscattering of seed NIR/UV photons causing the prompt 
optical flash. 



5.3 Dust extinction, luminosity, and energetics 

Figure [11] shows the optical afterglow curve of GRB 080603 A 
moved to a common z = 1 and corrected for the large 
dust extinction due to the si g htline within the host galaxy 
as described by iKann et al, ||2006|). The sam ple of other 
GRBs shown is taken from Kann et al.l l|2010l ). The after- 
glow of GRB 080603 A ranks among the mid/bright GRBs. 
Although not a dark GRB according to the /3ox < 0.5 defini- 
tion, this is a fair example of an optically observed dim burst 
mainly because of the large amount of dust within the host 
galaxy: the value Av,z = 0.8 mag is indeed, among those 
mea sured with good a ccuracy, one of the largest observed so 
far l|Kann et al.l|201 0'). This agre es with the findings from 
GRB host-galaxy studies (Perlev et al.l |2009| ) . samples of 



V 

GRBs with multi-colour photomet ric datasets (|Cenko et ahl 
|2009| ). and some individual GRBs (|Perlev et al.ll2011^ . Fur^ 
thermore, the SED we built is one of the very few which 
clearly favours an LMG2 extinction profile, with possi- 
ble evidence for the presence of the 2175 A bump that 
has rarely been observed in GRB afterglows ([Kriihler et al.l 
I2OO8I : lEUasdottir et al.ll2009l : iPerlev et ahlbom 

We could not directly measure the peak energy Ep of the 
prompt 7-ray spectrum; nevertheless, from the intermediate 
value of the photon index, F = 1.6, we can conservatively as- 
sume Ep^i = 160^130 keV. Combining this with the isotropic- 
equivalent radiated energy, Eiso = (2.2 ± 0.8 ) x 10^^ erg. 



GRB 080603A do es not violate the Ep-Eiso (|Amati et al.1 
l2002l : lAmatill201(il ) relation. 

Interpreting the late-time break as being due to a jet, 
we can provide an estimate of its opening angle, w hich in the 
ISM case turns out to be Oj = 5?7(-l?2, -fl?5) (|Sari et all 
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t (days after burst in the observer frame assuming z = 1) 

Figure 11. The optical afterglow light curve of GRB 080603A 
(thick line) moved to a common redshift of 1 and compared 
with the analogou s sample of long GRBs (grey lines) taken from 
iKann et al. I ll20ld) . All afterglow curves have been corrected for 
dust extinction. 

Il999h . We assumed standard values for the energy conver- 
sion efflciency, rj-y — 0.2 and for the particle density of the 
circumburst environment, n = 3 cm~^. The collimation- 
corrected released energy is Ej — (1.1 ± 0.4) x 10^° erg. 
This agrees well with the expect ations of the E^^i-E^ rela- 
tion (|Ghirlanda et al.ll2004l . |2007| ) , although the large uncer- 
tainty in i5p,i leaves this open. Calculating the analogous val- 
ues in the case of a wind profile (Nava ct al. 2006), these are 
very similar, although this density profile is disfavoured from 
the afterglow rise slope (Section lSX^ : Oj = 5?l(-0?9, +l°A) 
and E-y = (0.9 ± 0.3) x l(f° erg for the win d. Such an open- 
ing angle for the possi ble jet is v ery typical (|Zeh et al. I I2OO6I : 
iNavaet al...200& : .Racusin et al.|[2009r ). 



6 SUMMARY AND CONCLUSIONS 

GRB 080603A exhibits a number of properties which allow 
us to strongly test many aspects of the prompt and after- 
glow emission standard model. Our broadband dataset spans 
from the prompt 7-rays out to the radio band 13 days post 
burst, and also includes spectroscopic observations of the af- 
terglow as well as late-time photometry of the host galaxy. 
The main features of GRB 080603A are as follows: 

• a faint {R ~ 20 mag) optical fiash coincident with the 
last episode of a two-pulse, 150-s long 7-ray prompt burst; 

• a subsequent achromatic steep rise and peak around 
1600 s, which probably marks the afterglow onset; 

• no evidence for reverse-shock emission; 

• ISM circumburst environment favoured from afterglow 
modelling; 

• peak in the radio light curve detected at ~5 days, likely 
caused by the passage of the synchrotron spectrum peak; 



• late-time break in the afterglow light curve, interpreted 
as a jet break; the corresponding opening angle is 6j — 
5?7 (-l?2,-hl?5); 

• isotropic-equivalent 7-ray released energy _Eiso = (2.2 ± 
0.8) X 10^^ erg and a coUimation-corrected value of Ej — 
(1.1 ± 0.4) X 10^^° erg, both typical for long GRBs; 

• remarkable dust extinction within the host galaxy, 
A-v,z = 0.80±0.13 mag, that can be fit with an LMC2 profile 
(with marginal evidence for the 2175 A bump), and cannot 
be fit by the average MW and SMC curves, at variance with 
most GRB extinction profiles; 

• a comparable host-galaxy extinction (LMC; Av ~ 0.77 
mag) is required for fitting the host SED, possibly suggesting 
that the afterglow is being extinguished by a typical sightline 
through the host ISM; 

• extinction-corrected optical afterglow luminosity that 
lies in the mid-to-bright end of the distribution of GRBs at 
known redshift; and 

• projected offset from the host galaxy centre < 6 kpc, 
well within the offset distribution of long GRBs. 

Overall, the standard afterglow model seems to account 
for almost all of the observed properties of the broadband 
afterglow evolution. In particular, the best solution is given 
by an electron energy index distribution of p = 2, with both 
cooling and injection frequencies below the optical band at 
the time of the peak. However, the temporal evolution of the 
characteristic frequencies of the synchrotron spectrum can 
hardly be explained assuming typical values for the micro- 
physical parameters. 

We have constrained and crosschecked the Lorentz fac- 
tor in different ways: interpreting the optical afterglow peak 
as the fir eball deceleration y ields F = 130 ± 20. Secondly, 
following [Zou fc PiranI (|2010l ). we exploited the presence of 
a quiescent time between the two 7-ray pulses and derived 
an upper limit to F of 220. Finally, we focused on the optical 
and 7-ray prompt radiation to test whether inverse Comp- 
ton could be a viable mechanism to explain the GRB, and 
found that the allowed range for F is not compatible with 
the estimate deri ved from the afterglow properties (Fig. 1101 
IPiran et al.ll2009l ). Alternative interpretations of the optical 
flash, such as a reverse-shock origin, are excluded by the 
gap between the flash and the afterglow onset. Instead, an 
interesting possibility is that of an optical flash due to in- 
ternal shocks with a narrow distribution of ejecta Lorentz 
factors at larger radii. This would explain both the opti- 
cal flash being temporally disjoint f rom the afterglow onset 
and the peaky profile of the latter (|Meszaros fc Rees|[l999l : 
IPanaitescu fc Ve strand 201li ) . 

As an alternative to interpreting the optical afterglow 
peak as the fireball deceleration, the off-axis scenario re- 
quires an observer line of sight slightly off the jet cone, 
with (^obs ^ ^j) ~ 1°- The observed steep rise requires a 
structured jet with most of the energy in the jet core, with 
n ~ 4: for an angular profile modelled as £{9) oc {0/0j)~'' 
I Panaitescu fc Vestraiidll2008l ). 

Summing up, the clues derived for GRB 080603 A to un- 
derstand the nature of the prompt emission of GRBs high- 
light the importance of combining broadband campaigns 
and deep {R > 20 mag), rapid follow-up observations capa- 
ble of exploring the faint end of the prompt optical emission. 
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Table 8. Photometric data set of the NIR/optical afterglow of GRB 080603A. Uncertainties 
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^Midpoint time from the GRB onset time. 
''Corrected for airmass. 
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Table 9. Equivalent-width measurements for GRB 080603A 
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